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SUMMARY 

Three methods of measuring sound diffusion are described. One is derived 
from the curvature of sound decay curves and a second is based upon the variation of 
reverberation time with position. Both these yield useful indices of diffusion and are 
suitable for automatic analysis with the aid of a computer. A third method, using a 
rotating directional microphone, was found to be valueless for small rooms of moderate 
absorption, although it is reported to yield a useful index for large or very reverberant 
rooms. 



1. INTRODUCTION 

A sound field is said to be diffuse if: 

(a) there is uniform sound energy density within the 
region considered, and 

(b) there is equal mean energy flow in all directions at 
any point in the field. 

In the derivation of reverberation -time formulae, a 
diffuse sound field is assumed and sound absorption co- 
efficients must, strictly, be measured in a reverberation 
room in which the sound field is diffuse. However, current 
interest in diffusion measurement does not arise only from 
concern about the theoretical basis of reverberation time 
and the accurate measurement of absorption coefficients. 
The main reason is that certain acoustical defects in a room, 
e.g. flutter echoes, are associated with concentrations of 
sound energy in specific directions and that such defects 
can be cured by taking steps to scatter the energy more 
uniformly. 

Even now, reverberation-time and its variation with 
frequency are the only widely accepted objective criteria 
for the quality of a room; they are also widely recognized 
as necessary but not sufficient criteria. It has often been 
remarked that two rooms having the same reverberation- 
time/frequency characteristic can sound surprisingly dif- 
ferent; the determination of an index characterising the 
state of diffusion of the sound field in a room appears to be 
the next and, perhaps, the finat step in establishing satis- 
factory objective criteria. 

There are two difficulties about basing measurement 
techniques on the definition of diffusion. 

The first, a practical one, is that means of measuring 
sound energy and sound energy flow are not generally 
available. Microphones respond to sound pressure, sound 
pressure-gradient or a linear combination of these. While 
there are simple relations between sound pressure and energy 
density for the fields due to plane and spherical waves, in 
general no definite relations exist for the types of non- 
diffuse sound field found in rooms. The second difficulty, 
a theoretical one, arises from the fact that only the ideal 



state of perfect diffusion is defined quantitatively. 

Departures from the state of perfect diffusion can arise 
from inhomogeneity in the sound energy density and from 
anisotropic energy flow, so that a true measure of diffusion 
might involve two independent parameters. It is not sur- 
prising that most investigations into the measurement of 
sound diffusion have been based on empirical approaches, 
using parameters related only indirectly to the basic attri- 
butes of a diffuse field. Most of these parameters are 
unsatisfactory measures of diffusion in practice either 
because they are difficult to evaluate or because they are 
too insensitive. Others are applicable only to large rooms 
and are of limited value to a broadcasting authority having a 
large number of talks, interview and discussion studios. 

A practical study of sound diffusion in small rooms 
was made by Randall and Ward. They varied the surface 
distribution of absorbents and also used a number of dis- 
posable solid shapes which were attached to the walls and 
floor to scatter the sound on reflection. Parameters indi- 
cating the state of diffusion were derived by examining 
certain characteristics of the reverberation decay-curves. 
Useful though these parameters were, the time taken to 
evaluate them discouraged their general use. However, with 
the development of automatic reverberation time measure- 
ment, it was decided to investigate whether indices of dif- 
fusion based on the work of Randall and Ward could be 
obtained automatically from the digital decay data on 
punched paper tape. Two parameters appeared promising: — 

(i) P, a quantity representing the variation of slope 
of the decay curves with microphone position, 
was the standard deviation of the reverberation 
time, obtained from warble tone pulses, at 15 
microphone positions in the room. It was calcu- 
lated at seven frequencies from 500 Hz to 8 kHz 
and was found to vary markedly with changes in 
the distribution of absorbents. 

fii) S, a quantity related to the curvature of a decay 
curve, was the ratio of the slopes of the upper 
and lower halves of the decay curve. It was 
regarded as the most suitable criterion to be 
adopted for routine diffusion measurements in 
operational studios. 
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2. EXPERIMENTAL PLAN 

Randall and Ward showed that the state of diffusion 
could be controlled by a suitable distribution of absorbents 
over the surfaces of the room. Six different arrangements 
of absorbers were set up, in the experimental studio, which, 
from past experience, would be expected to cover the range 
of states of diffusion encountered in broadcasting practice. 
For each room condition, digital decay data on punched 
paper tape were obtained from the test recordings used in 
automatic reverberation time measurement. With suitable 
programming the paper tapes were then processed in a com- 
puter to yield parameters closely associated with P and S. 

In view of the complicated nature of diffusion it 
seemed likely that the two new computer-derived indices 
would not, at best, give precisely the same rank order for 
the six states of the room. At worst, they might not be 
capable of distinguishing the states of the room at all. As a 
check, it was decided to assess the states of the room by an 
independent method involving the use of a directional 
microphone. 

3. THE SIX STATES OF THE ROOM 

The experimental studio has the dimensions 6-7m x 
4-9m X 3-4m. The six different distributions of absorbents 
are shown in Fig. 1; for all conditions other than F there 
was a woolcord carpet with underfelt covering the floor. 
The absorbent surfaces consisted of a combination of low- 
frequency, mid-frequency and high-frequency absorbers in 
modular sections, of dimensions l-22m x 0-61 m, distributed 
evenly over each surface. 
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Fig. I - Distribution of absorbents for the six states of the 

experimental studio 

W/M Absorbent surface I I Reflecting surface 

The reverberation-time/frequency curve for each con- 
dition is shown in Fig. 2. Condition A would be acceptable 
for broadcasting purposes. 

The ten microphone positions for the reverberation- 
time measurements, and hence for the modified P and 5 
determinations, are shown in Fig. 3. The loudspeaker was a 
BBC Type LS5/5, pointed into a corner to reduce the level 
of the direct sound. An omnidirectional microphone was 
used. 

4. THE MEASUREMENTS 

4.1. Variation of Reverberation Time with Position 
The computer printouts of reverberation time for the 
six conditions enabled the standard deviation of reverbera 
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tion time at each frequency to be derived. When expressed 
as a percentage of the mean reverberation-time this gave the 
index P^, closely related to Randall and Ward's P. P^ was 
obtained with 1/3 octave-bandwidth noise-excitation from 
50 Hz to 8 kHz, at ten microphone positions, in contrast to 
P which was obtained with warble-tone excitation from 
500 Hz to 8 kHz at fifteen positions. 

The variations of P^ with frequency, for the different 
room conditions, are shown in Fig. 4. 

Since the variation of P^. below 1 kHz appeared essen- 
tially random, the mean P^ for the 1/3 octave-bands be- 
tween 1 kHz and 8 kHz was calculated for each condition 
and the results are shown in Table 1 . 



Room Condition 



4.2. Curvature of Decays 

First attempts to programme the computer to extract 
the index S from reverberation test-tape data produced 
absurd results. This arose from the difficulty of construc- 
ting an algorithm by which the 'upper half' and 'lower half' 
of a decay curve could be recognized. By dividing the 
dynamic range into two halves and programming the com- 
puter to run down the decay until the half-way level was 
exceeded, it was hoped that the mean reverbera tion -time 
above this point divided by the mean reverberation-time 
below, would yield the index S. Unfortunately, decay 
curves like that shown in Fig. 5 gave answers that were 
clearly wrong. 
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Fig. 4 ' Variation of reverberation time with position, P^. 
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Fig. 5 - Anomalous measurement of S 



Modifications to the definition of S were tried but 
they generally gave unreasonable answers with particular 
types of decay curve; that is, the indices derived were un- 
satisfactory because they were unduly sensitive to minor 
and acoustically unimportant characteristics of the decay 
curves. 

Fortunately, one modification to S was found to be 
satisfactory in this respect. The computer was programmed 
to construct a line passing through the commencement of 
the decay curve (point A in Fig. 6), with a slope equal to 
the mean slope of the decay down to 6 dB above the noise 
level (point B in Fig. 6). The ratio of the area bounded 
by the real decay curve AB and the lines AC and BC, to the 
area of the triangle ABC was taken as a measure of the 
curvature of the decay curve. The resulting index, called 
S^, was found to give results free from serious anomalies, 
presumably because minor characteristics of the curve make 
only a small contribution to the area under the curve. 

The data used for deriving S^ was the same as that 
used for P^. 
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Fig. 6 -Measurement of S^ 



It should be borne in mind that S differs from Sa in a 
number of ways: S was measured with pure tones and the 
results were averaged over twenty-six tone-bursts at 2 Hz 
intervals covering a band 25 Hz above and below the fre- 
quency of interest and at only one position in the centre of 
the room, whereas S^ was measured with noise of one- 
third-octave bandwidth at ten positions in the room. When- 
ever an individual reading of S exceeded unity the reciprocal 
of S was used in forming the mean of a set of readings, on 
the grounds that otherwise a set of convex and concave 
curvatures would tend to an average of unity. This pro- 
cedure was not followed in the determination of 5^ because 
it was felt that there is a greater danger that random errors 
in the measurement of S^ for a set of essentially straight 
decay curves would thereby be converted to systematic 
errors reducing the value of S^ to less than unity. 



The results of the measurements of S^ in the experi- 
mental studio are shown in Fig. 7. As with P^, there is 
little to distinguish the curves below 1 kHz. The mean 5"^ 
from 1 kHz to 8 kHz for each room condition is shown in 
Table 2. 
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4.3. Directional Microphone Measurements 

Meyer and others working at Gottingen, have develo- 
ped a method of measuring diffusion by examining the 
variation in the output of a highly directional microphone 
as it is rotated. While the method appears to have been 
successful for concert halls and other large rooms, it doer- 
not seem to have been applied to small rooms. 

Early attempts by the authors' colleagues to adapt the 
method for small rooms were fruitless,* apparently meaning- 
less results being obtained. Although a number of hypo- 
theses for the failure were put forward, to have tested them 
would have been uneconomic in the circumstances then 
prevailing. Instead, a more promising method was investi- 
gated. Furduev and Ch'eng T'ung have described a method 
whereby a directional microphone is rotated in a horizontal 
plane in a stationary noise field produced in the hall under 
consideration. The output of the microphone is plotted by 
means of a polar chart recorder, corresponding to the 
direction of the loudspeaker. The measurement is then 
repeated in an anechoic chamber, the input to the chart 
recorder being adjusted so that the response at 0° is equal to 
that obtained in the ball. 

In Fig. 8, R{6) is the polardiagram obtained in the hall 
and Did) is the directional characteristic obtained in the 
anechoic chamber; in the figure they are both surrounded 
by a circle of radius equal to the forward response of the 

• M.E.B. Moffat, L.M.B. Dyer and M.G. Croll devised most of the 
apparatus and carried out the early work. 



microphone. The index of diffusion d is defined as: 

A -A 
d--~l 



;i) 



where -4 is the area bounded by the circle and D(0) 
and A is the area bounded by the circle and R{Q). 




Fig.-S ' Furduev and Cheng T'ung differential-area polar- 
plot method 

Note that for a completely diffuse sound field R{0) 
would coincide with the circle and d would be unity. For 
a wave moving in one direction only, R[B) and D[6) would 
coincide making d zero. An attractive feature of Furduev 
and Ch'engT'ung's index is that it is amenable to theoretical 
analysis; they were able to deduce some interesting results 
regarding the indices obtained from microphones having 
figure-of-eight and cardioid polar diagrams. After studying 
the results they obtained in large rooms, using such micro- 
phones, it was decided to apply the method to the six 
states of the experimental studio: the object was to com- 
pare the results so obtained with the values of f^ and S^ 
obtained earlier. 

The rotating microphone turntable, an adapted gramo* 
phone turntable, is shown in Fig. 9. It was rotated at 
0-75 r.p.m. in synchronism with a B and K Type 2305 
tevel recorder fitted with a linear potentiometer. 

A typical chart record is shown in Fig. 10. On 
examining the charts for condition F of the experimental 
studio, it was decided that the method was impracticable 
in the form described by Furduev and Ch'eng T'ung. It 
was found that the point on the chart corresponding to 
zero input to the recorder was not at the geometrical centre 
of the chart but on a circle at about 1 0% of full scale from 
the centre. It proved valueless to measure the requisite 




Fig. 9 - Rotating microphone turntable in experimental 

studio 

areas on the chart with a planimeter, as recommended, since 
the effect of the offset zero could not be corrected for in a 
simple way, i.e. the diagrams would not scale simply. It 
was also evident that to have ignored the offset would have 
resulted in large errors. It would have- been possible to 
avoid this difficulty by plotting a linear chart rather than a 
circular chart, but values obtained in this way would not 
have corresponded to Furduev and Ch'eng T'ung's index. 

It was also difficult to ensure that the forward re- 
sponse always corresponded to a fixed radius. 

Improved results were obtained by using a digital form 
of microphone output using equipment previously develo- 
ped.^ The microphone output was amplified, rectified, 
sampled at V^ intervals, and the levels punched on eight- 
track paper-tape. The data was then processed in an 



Elliott 803B computer according to a programme for calcu- 
lating the requisite areas. 

Measurements were made with octave-bandwidth noise 
at 250 Hz and 500 Hz, and with 1/3 octave-bandwidth noise 
at 1 kHz, 2 kHz and 4 kHz. 

For each state of the room, two measurements were 
made at each of two microphone positions about 3-7 m 
from the loudspeaker. A PGS microphone (figure-of-eight) 
and a CI 2 microphone (cardioid) were both used. To get 
an overall picture, the values of d for the figure-of-eight 
and cardioid microphones at the frequencies 1 kHz, 2 kHz 
and 4 kHz were lumped together and the mean calculated 
for each condition; they are shown in Table 3. 



TABLES 
Room Condition A B C D E 

d 0-37 0-39 0-45 0-48 051 

5. SIGNIFICANCE OF THE RESULTS 

Fig. 11 shows a summary of the results by all three 
methods. The numbers in brackets next to the values of 
P^ and 5^ describe the rank order of the room conditions, 
assessed, on the basis of increasing P^ and decreasing S^. 



270 







Fig. 1 - Typical polar graph of the output of a cardioid microphone at 4 kHz in the experimental studio 



These two indices agree broadly in their assessments of the 
rank ordering of the states of diffusion of the room. The 
biggest discrepancy is in the rank order of condition C; 
experience suggests that P^, is the more accurate in this 
instance, as it ranks conditions B and D the least diffuse. 
These are the only two conditions in which one pair of 
opposite walls have absorption coefficients markedly dif- 
ferent from the other two pairs. 

It is difficult to assess the significance of differences 
in P^ because it is partly made up of errors in measurement 
of reverberation time. S^ does not suffer from this dis- 
advantage; for the experimental studio, the standard devia- 
tion of 5^, from 1 kHz to 8 kHz, ranged from 0-03 to 0-15. 
Since there were ten measurements of S^ at each fre- 
quency the standard error of the mean of S^ varies between 
about 001 and 0-05. 



The full value of P^ and .S"^ as diffusion indices can 
only be assessed after accumulating data and experience 
from a large number of rooms. Both indices have already 
been of some practical value, as was demonstrated when a 
reverberation test-recording from Studio 1 of the local- 
radio station at Leicester was processed. In the region 
around 400 Hz there was a marked increase in reverberation 
time, accompanied by a five-fold increase in standard devia- 
tion. The computer printout of Fig. 12 gives the corres- 
ponding values of 5^ , and shows marked curvature of decays 
in the same frequency region. Later investigation showed 
that a resonating metal cabinet was responsible. The very 
low values of S^ obtained suggests that poor diffusion 
caused by the presence of undamped resonating objects may 
affect 5^ much more than a poor distribution of absorbents. 
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Fig. 12 - Printout oj S^ for SiuJio I. Leicester 
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From Fig, 11 it is clear that the index d is of no 
value for small rooms of moderate absorption, since it 
increases with increase in reverberation time. It seems that 
in such rooms it is impossible to get far enough from the 
direct (and hence non-diffuse) sound from the source. The 
effect of this is simply to override any contributions from 
the directional characteristics of the reverberant sound field. 
The greater the absorption in the room, the smaller is the 
contribution from the reverberant sound. It is possible to 
argue that in general the more reverberant room is ipso facto 
more diffuse, but diffusion indices based on this view are 
not of very great value because they do not correspond to 
subjective judgements of the effects due to the poor distri- 
bution of absorbents. 

An attempt was made to separate the contribution of 
the reverberant field by assuming a relation of the form: 



rf^b^+b^.r+b^j', +b^5^ 



Fig. 11 - Comparison of the three indices of diffusion 



where the b's are constants and T is the reverberation time. 
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The sets of values of d, T, F^ and S^ were run in the 
Elliott 803B computer with a programme designed to deter- 
mine the influence of each of the variables.* The value of 
b| obtained was significantly positive, but the values of b 
and b were not significantly different from zero. The 
attempt therefore failed, but of course the result does not 
exclude the possibility of relations other than the linear one 
postulated above. 



A diffusion index due to Furduev and Ch'eng T'ung, 
which is based on the measurement of areas of polar dia- 
grams obtained with a directional microphone, is found to 
be of no value in small rooms of moderate absorption. 



6. CONCLUSIONS AND RECOMMENDATIONS 

Modifications have been made to the diffusion indices of 
Randall and Ward, P and 5, to make them suitable for 
derivation from computer analysis of reverberation test- 
recordings. 

The new parameters, P^ and 5^, like the old, depend 
on the variation of reverberation time with position and the 
curvature of decay curves, respectively. P^ and S^ also 
have been shown to be dependent on the state of distribu- 
tion of absorbents. They can be obtained very easily from 
the computer analysis of reverberation test-recordings with 
little extra expenditure of time and effort. Their full sig- 
nificance will only be realized when they have been deter- 
mined for a sufficient number of operational studios. It is 
recommended, therefore, that P(. and S^ be determined and 
noted whenever reverberation test-recordings are made. 



Library programme S.101, 'Multiple Regression'. 



7. REFERENCES 

1. RANDALL, K.E.a^(i WARD, F.L. Diffusion of sound 
in small rooms. Proc. I.E.E., 107B, 35, pp. 439 - 450, 
(Sept. I960}. 



2. Automatic reverberation time measurement. 
Research Department Report in preparation. 



BBC 



3. RICHARDSON, E.G. fl«i/ MEYER, E. (Eds.) 
aspects of sound —Vol. 111. Elsevier, 1962. 



Technical 



4. FURDUEV, V.V. and CH'ENG T'UNG. Measurement 
of diffuseness of the acoustic field in rooms by the 
directional microphone method. Soviet Physics — 
Acoustics, 6, 1, pp.' 103 - in, July -September 1960. 

5. An analogue to digital conversion system and paper tape 
transcoder for processing acoustic signals, BBC Research 
Department Report No. PH-8, Serial No. 1967/29. 



SMW 



Printed by BBC Reseorch Department, Kingswood Warren, Tadworlh, Surrey 



